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Starting from the inhomogeneous shear-free Nariai metric we show, by solving the Einstein- 
Klein-Gordon field equations, how a self-interacting scalar field plus a material fluid, a variable 
cosmological term and a heat flux can drive the universe to its currently observed state of homo- 
geneous accelerated expansion. A quintessence scenario where power-law inflation takes place for a 
0^ ' string-motivated potential in the late-time dominated fleld regime is proposed. 

a^ 

I. INTRODUCTION 

The theory of general relativity implies that the large scale geometry and evolution of the Universe is dictated by 
,_^ • its content of matter and fields, including energy fluxes, shear stresses, particle production, and so on. Recently, there 
\ have been claims in the literature that the Universe, besides its content in normal matter and radiation, must possess 
C\| ■ a not yet identified component (usually called quintessence matter, Q-matter for short) ||l|, [gj, |g], Q, characterized 
^ ' by a negative pressure, and possibly a cosmological term which may be constant or not [pi, [pi. These claims were 
^^ prompted at the realization that the clustered matter component can account at most for one third of the critical 
^_J density. Therefore, an additional "soft" (i.e. non-clustered) component is needed if the critical density predicted by 

vQ many inflationary models is to be achieved. 

f^ ' Very often the geometry of the proposed models is very simple, just Friedmann-Lemaitre-Robertson- Walker 

^\ , (FLRW), partly because of mathematical simplicity, and partly because they deal with not very early stages of 

On ' cosmic evolution. However, up to now no convincing argument has been avanced to the effect that the geometry of 

^j ' the very early Universe (say, some time before nucleosynthesis) could not have been either anisotropic or inhonio- 

^^. geneous, or both. Note that it is very natural to assume the geometry at that primeval epoch more general than 

J ' just FLRW. Moreover, recently it has been demonstrated that given any spherically symmetric geometry and any 

(2J[)i set of observations, evolution functions for the sources can be found that will make the model compatible in general 

IV with observation [Q . Further motivations to study the evolution of inhomogeneous cosmological models can be found 

in references [g| and [Q|. In constrast to FLRW models, inhomogeneous spaces are in general compatible with heat 

fluxes, and these might imply important consequences such as inflation [nffl or the avoidance of the initial singularity 

V^ [0| . Here we focus on an isotropic but inhomogeneous spherically symmetric universe which besides a material fluid 

contains a self-interacting scalar field (which can be interpreted as Q-matter), and a cosmological term, A which, in 

general, may vary with time. As it turns out the homogeneization of the Universe at large times depends very much 

on the value the effective adiabatic index 7, defined below (Eqn. (|5|)), takes. When more generally this quantity 

varies with time, homogeneization can also be achieved under rather ample conditions. 

Obviously, density inhomogenities triggered by gravitational instability must be present at any stage of evolution, 
even when the universe has reached a near homogeneous state describable in the mean by the FLRW metric -i.e. an 
average homogeneous universe housing growing as well as decaying density perturbation modes, the former leading 
eventually to the formation of the cosmic structures we observe today. About the evolution of scalar (i.e. density) 
inhomogeneities there is a rather large and still growing body of literature (see e.g. |12| and references therein) but we 
shall not deal with them as they lie outside the main focus of this paper. We only mention that the negative pressure 
associated to Q-matter and A will tend to slow down the growing modes (see e.g. [y, |1^, [|l^), and shift the epoch 
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of matter-radiation equality toward more recent times jlj]. A detailed study of all this will be the subject of a future 
work. 

Section II presents the basic equations. Section III solves the Einstein-Klein-Gordon (EKG) field equations assuming 
7 (but not A) a constant. Section IV studies the asymptotic evolution toward a Q-matter dominated era. Section V 
solves the EKG equations when 7 is a function of time. Finally, section VI summarizes the findings of this paper. 
Units have been chosen so that c = G = 1. 

II. EINSTEIN-KLEIN-GORDON FIELD EQUATIONS 

Let us consider a shear-free spherically-symmetric spacetime with metric p5[ 

ds^ = -A{t, rf dt^ + B{t, rf [dr^ +r^ dfl^] , 

where as usual dfl^ = d6^ + sin^ 9 dcjp' . For later purposes we introduce the functions F{t,r) = 1/B{t,r) and 
v{t,r) = AF. Thus, the line element assumes the more convenient form 

ds^ = ^— ^ \-v(t,rf dt^ + dr^ + r^ dn^] . (1) 

F(t^ ry 

As sources of the gravitational field we take: a fluid of material energy density pf = pf{r,t), hydrostatic pressure 
Pf ~ P/(''jO; with a radial heat flow (q^ = qr{r,t) and qt = qe = q<f, ^ Q), plus a cosmological term, related to the 
energy density of vacuum by A = ^Trpyac, that depends only on time A = A(i), and a self-interacting scalar field 
(j) driven by the potential V{(j)). Taking the scalar field to depend only on t, its energy-momentum tensor may be 
written in the perfect fluid form 

Tik = (p0 + P4,)uiUk + P4,gik , (2) 



where m* = (t)'"^ / y/—(j).i4' together with 

P0 = -^'/',i0''-n0)- (3) 

The fluid interpretation of the scalar field, while not compulsory, has proven very useful in the study of the inflationary 
phase and reheating of the universe -see for instance M& and p7| . In particular it leads to consider its equation of 
state Pel, — (70 — 1)P0- Hence the scalar field can be interpreted as Q-matter since, depending on V{(j)), 7^ may be 
lower than one -see e.g. [pi. The stress energy-tensor of the normal matter, with a heat flow, plus Q-mattcr (scalar 
field) and the cosmological term is 

Tl = (Pf +P4>+Pf+ P<pWuk + (A - Pf - p^)dl + q'uk + qku' , (4) 

where the comoving four-velocity is normalized so that u^Ui = —1, with u* = 5l{gtt)~^''^, and hence u* = F/v. 
Obviously the heat flow is orthogonal to Ui (i.e. q^Ui = 0). As equation of state for the fluid we choose p/ — (7/ — 1) pf 
where 7/ is a function of t and r. Taking into account the additivity of the stress-energy tensor it makes sense to 
consider an effective perfect fluid description with equation of state p = (7 — 1) p where p = pf -\- p^^ p = pf + p^ and 

^^IfPl + l^^ (5) 

Pf + P4> 

is the overall (i.e. effective) adiabatic index. 

The requirement that the cosmological term A is just a function of t leads to the restriction that 7 also depends 
only on t to render the system of Einstein equations integrable. The nice result we are seeking is a solution that has 
an asymptotic FLRW stage, with A evolving towards a constant, and the heat flow vanishing in that limit. 

For the metric (H) the EKG equations are 

p2 
12xF' -12 FF' -8 xFF" -3^^ + p + A = 0, (6) 

V 



p2 FFv FF F^v' 
12xF'^-8FF'-3^~2^-^ + 2^+A^-^ 



y2 y3 y2 



-p + A = 0, (7) 
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,2 



F^ „FFv 
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FF xFF'v' F^v' xF^v" 
2^p-8 +4 +4 p + A = 0, (8) 



v'^ V V V 



^xF^Fv' ^xF^F' 
- 8 5 + 8 + 9r == 0, (9) 

V V 

0-2f^+i,f = O, (10) 

where x = r^, an over-dot means d/dt, and a prime d/dx. Wc thus have five equations and seven unknowns 
{F ,v , p ,p, (j), Qr and A). Hence two more equations are needed to render the system determinate. These are the 
above equation of state relating p and p and the condition dA/dx = 0. 
From (0) and (||) we get 

v" F" , , 

then using the ansatz F = [a{t) + b{t)x] and v = [c{t) + d(t)x]^ in ( [ll| ) we find that the constraints n (n — 1) = 
2k {k — 1) and a{t) d{t) = b{t) c{t) must be satisfied. Next we study the simplest case, namely n = k = 1. This is the 
most interesting instance because the functions a{t),b{t),c{t) and d(t) can be chosen freely. Another set of solutions 
can be obtained if the relationship Vxx/v = 2[i^(i)]^ is assumed in (|l l[) . These solutions are presently under study and 
we shall report on them in a future paper. 

A. Einstein equations with time-dependent A and 7 

When n = A: = 1 we obtain a set of solutions that contains those of Modak (6 = 0) Il8], Bergmann (c = a, d = &) 
II19I and Maiti (b — d — k aJA , with fc = 0, ±1) [g^- Another possibility arises when d ^ 0. This solution can be also 
obtained from the metric p) by imposing this metric be conformal to Minkowski's. Namely, the integrability of the 
transformation —vdt + dr = drj and vdt + dr = da leads to u = v(t) and n ^ k = 1. This allows us to re-define the 
time by vdt -^ dt. Then the metric (|l|) becomes 

ds^ = ^ ^ (- df + dr^ + r^ AV?\ , (12) 

and equations (|6|), ([7|) and (|9|) turn into 

p + A= 12a6 + 3a^ + 6a5a; + 3fe^a;2, (13) 

p- A= hb'i) - 36^) x^ + 2 hb"^ - 3db + ab + db\ x - Sab - 3d'^ + 2ad , (14) 

and 

qr = -4,^/xb{a + bxf , (15) 



respectively. We next impose that A depends solely on time, so it must have the form 

A{t) = 12ab + 3a'^ + f{t), (16) 



where f{t) is a function to be determined. Then, (|13|) and (14) imply 



p{x,t) =6dbx + 3b^x^ ~ f{t), (17) 



p{x, t) = (266 - 36^ j x^ + 2 [2b^ - 3ab + ab+'dbjx + iab + 2ad + fit). (18) 

We further impose the equation of state p — (7 — 1) p. Taking into account that a and 6 depend only on time, and 
equating the coefficients of same power of x we obtain a set of equations to determine a, 6 and / 

bb-^jb' = 0, (19) 

«-37ja+Y^a = -26, (20) 

/ = --(26 + a). (21) 

7 

To show the conditions that lead asymptotically to a FLRW metric it is expedient to introduce the time coordinate 
dr — dt/a. Thus dl^) becomes 

ds^ = ^- 2 [-dr^ + i?2 (dr^ + r^ dn^)] , (22) 

(1 + Mr^) 

where M = b/a and R = l/|a|. This metric is conformal to FLRW, and the conformal factor approaches unity when 
M^O. 

III. CONSTANT ADIABATIC INDEX 

When 7 is a constant different from 2/3, the general solution of ( [l9|) and (|2C| ) becomes 

b{t) = KAt^ , (23) 

a(i) = CiAt'^T^ + CsAt"^ - -KAt^^ , (24) 



thereby 



Whereas for 7 ~ 2/3, we get 



and 



M(i) = -ft: ( Ci + CaAi^i - -KAtA . (25) 



6(i) = Ke^^\ (26) 

a(i) = (Ci + CaAi - if At^) e^'^*, (27) 

M{t)=K{Ci+C2At-KAt'^y\ (28) 



where K, C, Ci and C2 are arbitrary integration constants, and /Si — t — tQ with io some initial time. Inserting (|2J), 
(H), (ll) and @ in (|l|), dll) and (|l|) we get 






and 



qr{r,t) 



SKrAt"^"^ r 



9 (37 - 2) 



for 7 7^ 2/3, and 



3C2Ai"^5^7^ + 3 (Ci + ifr^) At'W^ - KAt^^^^ 
A{t) = (3C2^ + 12XCi) e^c"-^* , 



q,{r, t) = -4if r (Ci + CaAf + ^(r^ - At^))^ e^^^* , 



(29) 



(30) 



(31) 



(32) 



for 7 = 2/3. In the following we analyze the asymptotic behavior of these solutions in the limits for t such that 
M -^ 0. We note by passing that in the M — > limit the time coordinate r becomes the cosmological time, as can 
be seen from (p2). Two alternatives of asymptotically expanding universes appear depending on the map between t 
and T -equations (||), (||) and (||), (||) below. 

Case At -> 



In this limit we obtain 



CzAt^^^ , 



(33) 



2 — 37 A g '^-^T 

At ~ — -Ar^^^ 

2C2 (1 - 37) 



(7.^1/3), 



(34) 



At~— iTiAt (7 = 1/3), 

L^2 



(35) 



R{t) 



C2 



Rir) 



2C2 (1 - 37) 



2-37 



At 



(7 ^ 1/3), 



-CaAr 



Co 



(7=1/3), 



(36) 

(37) 



A(r) 



3(2-37)^ 

4(l-37)^Ar2 



(7 7^ 1/3), 



A ~ 3Cl (7 = 1/3), 



(38) 
(39) 



, , SKC2 
g.(r,r)^^— ^r 



2C2 (1 - 37) 
2-37 



At 



(7 7^ 1/3), 



(40) 



qr{r,T) -^SKC^re 



2^„3C2At 



(7 = 1/3). 



(41) 



When 1/3 < 7 < 2/3 we have, for large cosmological time t, an accelerating universe that homogenizes with vanishing 
cosmological term and heat flow. In this stage we may define an asymptotic adiabatic index by equating the the 



exponent in equation ( |3q ) to 2/(37asy), i.e. 7 = 4/3 (4 — 37asy)- In this range one has 1 < 7asy < 00 leading to a final 
power-law expansion era. For 7 = 1/3 the map between t and r changes and we have an asymptotically de Sitter 
universe with finite limit cosmological term. For the remaining values of 7 the universe begins at a homogeneous 
singularity with a divergent cosmological term. When 7 < 1/3, the heat flux asymptotically vanishes near the 
singularity, while for 7 > 2/3 it diverges. 

Case At ^ 00 

In this limit we obtain 



1 c "y — ^ 

3 



(42) 



Ar 



3 2 



if 4 -37 



37 . -4.+31 



(43) 



Rir) ^ -^ 



K{A- 37) 
3(2-37) 



Ar 



(44) 



A(r) 



24(2-37)^ 
(4 - 37)^ At2 '' 



(45) 



g,(r,r)~ -24 



for 7 7^ 2/3, and 



(2 - 37)^ r 
(4-37)^ A73' 



a^-i^At^c^^*, 



At =:. e-^^*, 



(46) 

(47) 
(48) 



R{t) ~ -CAr, 



(49) 



A(r) 



(3C| + 12ifCi) 

A^ 



(50) 



9r(r,T) 



7«T-^ (InAr)^ 



(51) 



for 7 = 2/3. 

In this case the asymptotic adiabatic index is given by 7 = (97asy — 4)/(97asy — 3) and we have two alternatives 
for asymptotically expanding universes. When 2/3 < 7 < 1 the universe homogenizes for large cosmological time 
with vanishing cosmological term and heat flow. We note that even though an asymptotic negative cosmological term 
occurs, the universe ends in a power-law evolution R oc Ar" with < a < 1. When 7 = 1, the late time evolution 
changes to an asymptotically Minkowski stage. For 1 < 7 < 4/3 the universe starts homogeneously in the remote 
past with a vanishing scale factor, cosmological term and heat flow. For the remaining values of 7 the universe begins 
at a homogeneous singularity with a divergent cosmological term. 

An exact solution with explicit dependence on the asymptotic cosmological time r can be found when the integration 
constants Ci and C2 vanish. In such a case each approximate expression (H2h-(B6|) becomes an equality and the metric 
reduces to 



ds' ^ 



1 



2 2-37 



-dr^ +At^^^ ( di 



'dO") 



(52) 



where m is a redefinition of the old integration constant K, the adiabatic index 7 and tq. The last constant was 
introduced by scaling the radial coordinate r —^ ror. 

Proposed varieties of soft matter with 7 < 1 include cosmic string networks |21|], "K~matter" [g2[ and quantum 
zero-point field [g3|, all with 7 — 2/3, as well as "QCDM" (for "unknown cold dark matter") with ~ 0.4 [Q and the 
quintessence (or Q-component) with < 7 < 1 B. Fluids with values of 7 less than 2/3 may be termed "inflationary 



matter". Equation (36) shows an evolution that corresponds to this kind of matter. 

These results illustrates how homogeneization of a universe dominated by matter that has negative pressure in the 
present era may have occurred. A smooth unclustered dark matter component with negative presure could reconcile 
a flat, or nearly flat, universe with a density in clustered matter well below the critical value, and moreover explain 
the recent high redshift supernovae data suggesting that the universe is currently under an accelerated expansion [g4| , 
p5[ . For a perfect fluid negative presure leads to instabilities that are most severe on the shortest scales. However, 
if instead the dark matter is a solid, with an elastic resistance to pure shear deformations, an equation of state with 
negative presure can avoid these short wavelength instabilities. Such a solid may arise as the result of different kinds 
of microphysics. Two possible candidates for a solid dark matter component are a frustrated network of non-Abelian 
cosmic strings or a frustrated network of domain walls. If these networks settle down to an equilibrium configuration 
that gets carried along and stretched by the Hubble flow, equations of state result with 7 = 1/3 and 7 = 2/3, 
respectively. One expects the sound speeds for the solid dark matter component to comprise an appreciable fraction 
of the light speed. Therefore, the solid dark matter does not cluster, except on the very largest scales, accessible only 
through observing the GBR anisotropy at large angles |j]. 

IV. ASYMPTOTIC EVOLUTION TO A QUINTESSENCE-DOMINATED ERA 

As a first stage of towards more general scenarios with a slowly time-varying 7, we will explore a model that evolves 
towards an asymptotic FLRW regime dominated by Q-matter (i.e. the scalar field). We will show that this system 
approaches to the constant 7 solutions for large times found above. In this regime equations (g) and (nw become 

3ij2~p^ + 1^2^ ]/(</,)+ A, (53) 

+ 3i?^+^^O, (54) 

where H — R/R and a dot means d/dr in this section. /^From these equations and (|^) it follows 

and 

7* = T--^ • (56) 

* ^02 + n0) 

In last section we found that the general asymptotic solution for the scale factor R{t) oc At" has the power-law 



behaviors (36), (|44|) for any value of the effective adiabatic index 7. Then, using these expressions and (|38|) and (|4E 



together with (B3[) and (pa) in (p6[), the latter becomes 



where 






(58) 
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3(2-37)' 

P = 91 

4(1-37)' 


.1 2, 

^3<^<3)' 




2(37-1)' 


a = 


6(1-7) 
4-37 ' 


24(2-37)' 
(4-37)' 


(|<7<1). 



(59) 



To investigate the asymptotic limit in which the energy of the scalar field dominates over the contribution of the 
perfect fluid we assume that 3ajf > 2. In this regime the two terms in (|57| ) proportional to the energy density pf are 
positive and neglegible. The adiabatic scalar field index can be approximated by 



^^"3^ 



1 



2 



(60) 



where a — Pf/p^ <S^ 1. Inserting these equations in (||) we obtain the first correction to the effective adiabatic index 
for the solutions (p8|) 



1± J(37/-2)a 



(61) 



and 7 ~ 2/3 + 7/cr/2 for the solution (|59|). The negative branch of ( |6l| ) yields a consistent asymptotic solution for 
the range ^ < 7 < |. We note that this solution describes a deflationary stage with a limiting exponent a = 1. 

Oftenly power-law evolution of the scale factor is associated with logarithmic dependence of the scalar field on 
proper time |2q ]. Thus, assuming that (f>{T) ~ Clnr with the constant C to be determined by the system of equations 
(p3|) and (Isl), and using these expressions together with (|38|) and (|4|) in (|53|) and (54) it follows 



3a^ 

^2 



2t2 



V 



(3 



(62) 



(3a - 1) C dV ^ 
From (p^), ( |63| ) we obtain the leading term of V{(1>) for large (/; 



and 



for i < 7 < 



A' 



37, 



Vo 



|, while 



A' 
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17-427 + 2772 ' 



Vo 



37 + 2 
37(37-1)' 



2 (-14 + 15 7) (17 - 42 7 + 2772 
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(63) 



(64) 



(65) 



(66) 



for I < 7 < 1. Inserting the dominant value of the effective adiabatic index in (pa), (pw we find ^2 = 2, Vq = 2 and 
C=l/V2. 

The models considered in this section are based on the notion of "late time dominating field" (LTDF), a form 
of quintessence in which the field (j) rolls down a potential V{(j)) according to an attractor solution to the equations 
of motion. This solution is an attractor since for a very wide range of initial values for M, (p a-nd (j) it rapidly 
approaches a common evolutionary path, i.e. the late behavior is insensitive to the initial conditions. This model 
has an advantage similar to inflation in that for a wide range of initial conditions the universe is driven to the same 
final evolution. The ratio a of the background fluid to the field energy changes steadily as (j) proceeds down its path. 
This is desirable because in that way the Q-matter ultimately dominates the energy density and drives the universe 
toward an accelerated expansion. 

Recently Ferreira and Joyce P7| proposed a model based on an exponential potential. Their self-adjusting solutions 
are attractors and U,^ remains constant for a constant background equation of state as 7^ = 7/. {^4, changes slightly 
when the universe shifts from radiation-dominated to matter-dominated expansion). This means, for example, 
that Vi^ is constant throughout the matter-dominated epoch. For a constant Q.^ to satisfy the structure formation 
constraints requires Vi^ < 0.2 and Vtf > 0.8. This however runs into conflict with the current best estimates of Vlj 
and produces a decelerating universe at variance with recent supernovae observations 124], | p5| . 

By contrast, our LTDF model only requires that the potential has an asymptotic exponential shape for large (j). 
So, the interesting and significant features of our model are: (a) like the self-adjusting case, a wide range of initial 
conditions are drawn towards a common evolution; however, (b) the LTDF solutions do not "self-adjust" to the 



background equation of state but rather, maintain some finite difference in the equation of state such that the field 
energy eventually dominates and the universe enters a period of acceleration. Compared to the self-adjusting model, 
ours does not require any additional parameter and allows a much wider range of potentials, provided they have an 
exponential tail. LTDF solutions exist for a very wide class of potentials. The energy density of the field decreases 
as i?~^, and this power-law remains unchanged in every epoch of the universe when the background expansion turns 
from radiation- to matter- to quintessence-dominated. The value of 7^ differs from the background equation of state 
such that the value of fJ^ increases as the universe ages. Hence, i7^ grows to order unity late with time. 



V. GENERAL SOLUTION FOR A TIME-DEPENDENT ADIABATIC INDEX 



A time-varying 7 is very natural because different matter components redshift at various rates. Then the question 
arises whether the conditions leading to homogeneization we have found with constant 7 also hold when 7 varies with 
time. The rationale behind this approach is the following. When different components enter the stress-energy tensor 
of the cosmic medium, it is natural to expect that each epoch is dominated by the energy density of just a single 
component with a constant, or nearly constant, adiabatic index, say 701. The others components can be regarded as 
small perturbations. As the universe expands, sooner or later, the component with the adiabatic index inmediately 
lower than 701 (say 702) takes over, and a new epoch begins. 

We start by giving the general solution to equations (|l^) and (gO|) when 7 is a function of time 



exp 



dtw 



(67) 



where 



a = —2 exp dtw] dtw 



dt 



(68) 



/di(2-37)' 
provided 7 ^ 2/3. Inserting (|8|) and (|6^) in (|l|), (H) and (|l^) it follows 



(69) 



A = AK'^ exp 2 dtw 



w'' J 7 



dtw-^ I i,] /4 



2w^ 

1 



dtw'^ I ^r] - 6 / di wM ^r 
W'' / I I w^ 



(70) 



qr{r,t) — —AKrwexj) 3 / dtw 



-2 I dtw'' I -^+r' 



(71) 



Though we have reduced the general coupled system of equations to quadratures, it is very involved to obtain the 
solution in closed form except for constant 7. This is why next subsection focus on approximated solutions assuming 
that 7(t) admits an analytic expansion. 



A. Homogeneization with a varying adiabatic index 



The transition previously described from one epoch to the next, assumed to be gentle, may be modeled by an 
homographic function 7(i) = {At + B)/{Ct + D) where A, S, C, andD are constants [|8|. With this in mind, we next 
investigate two kind of behaviors for the fluid components that lead to a final homogeneous stage characterized by 
M{t) — > 0. These can be associated with two asymptotic evolutions of 7(t). 



1. In the limit t ^ we assume that 7(t) has a Taylor expansion 

with 701 = B/D and 711 = {AD-'^ - BD-"^). Then, for 701 ^ 2/3 wc find that 



711 



A 
D 



B 



M{t) = 



Kt 

'c7 



1 



2 - 3701 



\nt + 0{t^) 



and that the first corrections to equations ( p3|) and ( [40|) vanish as tint, while the first correction to 
higher order. 



2. In the limit i — *■ 00 we assume for j{t) the expansion 



l{t) = 702 + ^ + O{t-^) 



702 = A/C , 712 = ^ - ^ , 



obtaining for 702 7^ 2/3 



M{t) = - 



f2 



67 



12 



Ini 



2-3702 t 



+ o(t-i) 



(72) 

(73) 
|) is of 

(74) 
(75) 



In this case the first corrections to equations ([42|), (45) and ([46|) vanish as Int/t 



These results show that homogeneization occurs under the same conditions for both 70 as it does for constant 7 for 
a wide range of evolutions of the adiabatic index, provided it has a constant limit, 701 and 702 for t —f and t —>■ 00, 
repectively, and is analytic about these points. 

VI. CONCLUDING REMARKS 

We have investigated a class of solutions of the Einstein field equations with a variable cosmological term, heat 
flow and a fluid with variable adiabatic index that includes those of Modak, Bergmann and Maiti and contains a 
new exact conformally flat solution. We have also considered the contribution of a homogeneous minimally coupled 
scalar field to the stress-energy tensor. The solutions of the EKG system that lead to an asymptotic FLRW stage 
were analyzed when the adiabatic index remains constant. We have found that asymptotically expanding universes 
occur when 1/3 < 7 < 1 that homogenizes for large cosmological time with vanishing cosmological term and heat 
fiow. For 1/3 < 7 < 2/3 the evolution is given by ( p6| ) and corresponds to a power-law accelerated expansion for 
large cosmological time r, as follows from (^4|). On the other hand, when 2/3 < 7 < 1 even though an asymptotic 
negative cosmological term occurs, the universe evolves toward a decelerated expansion. The particular case 7 = 1/3 
leads asymptotically to a de Sitter universe with a finite limit for A. We have shown that these results also apply to 
a time dependent adiabatic index. 

We have carried out a detailed analysis of a model in which Q-matter dominates over cold dark matter. This LTDF 
solution is an attractor because, even for large initial inhomogeneities and a wide range of initial values for cj) and (j), 
the evolution approaches a common path. It was shown that this model can be realized for a wide range of potentials 
provided they have an exponential tail. This is quite interesting because recently, there has been increasing activity 
related to scalar fields with the Liouville form (exponential) potential. It arises as an effective potential in many 
theories such as Jordan-Brans-Dicke theory, Salam-Sezgin theories and superstring theories. Indeed, most theories 
undergoing dimensional reduction to an effective four-dimensional theory, result in a linear combination of exponential 
potentials, and one of these will eventually dominate. 

This work can be generalized in different directions. The more obvious one is to allow A and 7 to depend on position 
too, but this seems analytically impracticable. A less hard possibility is to consider non-flat spatial sections. This is 
interesting since although the location of the flrst acoustic peak in the angular power spectrum of the GBR suggests 
a flat universe, its exact position is still uncertain [ p9[ . Therefore it may well happen that the universe is open or 
even closed. Again the difficulty with the corresponding analysis is basically mathematic since it is doubtful that in 
this more general case the EKG field equations admit analytical solutions. A seemingly less involved generalization is 
to include a bulk dissipative pressure in the stress-energy tensor ^ , something rather natural [Q ■ Its effect should 
tend to further accelerate the expansion, as it lowers the total pressure of the cosmic fluid. The problem with this is 
the lack of fully realiable expressions for the coefficient of bulk viscosity of the fluids involved in the hydrodynamic 
description. However, in this regard the situation is no worse than that encountered in inflationary scenarios in which 
potentials for the scalar field are frequently proposed with no obvious physical ground. 
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